On &;-Core Percolation in Four Dimensions 
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The fc-core percolation on the Bethe lattice has been proposed as a simple model of the jamming 
transition because of its hybrid first-order/second-order nature. We investigate numerically fc-core 
percolation on the four-dimensional regular lattice. For k = 4 the presence of a discontinuous 
transition is clearly established but its nature is strictly first order. In particular, the fc-core density 
displays no singular behavior before the jump and its correlation length remains finite. For fc = 3 
the transition is continuous. 



After its introduction fc-core percolation has been 
proposed to be relevant in a variety of different contexts, 
see 0, H| • The problem, also referred to as " Bootstrap 
percolation", is defined as follows. The sites of a given 
lattice are populated with probability p. Then each site 
with less than fc neighbors is removed, the procedure be- 
ing iterated until each site has at least fc neighbors. 

On regular lattices in dimension d the model exhibits 
two different behaviors depending on the value of fc. If 
k > d the cluster must be extended in order to survive 
the culling process but it is completely decimated for any 
p < 1 in the large size limit £| . The behavior for large but 
finite system size has also been investigated 0,0 Hill but 
strong disagreement between the theoretical predictions 
and numerical simulations has been found. The highly 
non trivial origin of this discrepancy has been clarified 
only recently |j ■ If k < d it is easy to realize that there 
are small "self-sustained" structures (e.g. d-dimensional 
hypercubes) that can survive culling irrespective of their 
environment. In this case the fc-core always exists and 
the problem is rather if it percolates or not and what is 
the nature of the percolation transition [ljj . 

On the Bethe lattice the fc-core percolation transition 
is known to be discontinuous Starting from high 

values of p the density of the fc-core drops discontinu- 
ously to zero at p c . The transition however is not simply 
first order, the density near the transition is given by 
p(p) ~ p(p c ) + b(p — pc) 1 / 2 . Furthermore it has been re- 
cently pointed out that the transition is also accompanied 
by diverging correlation lengths. This behavior has mo- 
tivated the proposal of fc-core percolation as a model of 
the jamming transition ^^l- There are indeed evidences 
that this transition has a mixed character 0] . 

This brought new interest into the question of wether 
the hybrid nature of the transition in the Bethe lattice 
survives in finite dimensions. This is certainly not the 
case for cubic lattices in d < 3. Indeed for fc = 2 the 
transition is continuous and has the same critical point of 
ordinary percolati on fT3 . For fc = 3 and d = 3 the transi- 
tion is continuous |lllll5l Il6t ITt| with exponents consis- 
tent with those of ordinary d = 3 percolation [l6j . These 
results are valid on cubic lattices and they not exclude 
the possibility of a mixed transition for d < 3 provided 
the structure of the lattice or the constraints are differ- 
ent. Indeed recently a 2-dimensional model with a mixed 



transition was exhibited 0] and numerical evidences of 
a mixed transition in another 2-dimensional model were 
reported in |l2j . As for regular lattices, an expansion 
in powers of 1/d has proven that turning on dimension 
perturbatiyelv does not destroy the mixed nature of the 
transition [l(| , thus suggesting that the hybrid transition 
may exist for some (c?>3,2<fc<cJ+l). In this work 
we investigate numerically fc-core percolation on the four- 
dimensional hypercubic lattice. For fc = 3 we have found 
a continuous transition and we did not further investigate 
the critical behavior. In the case (d = 4, k = 4) we find 
negative results concerning the hybrid transition: while 
the presence of a discontinuous transition is clearly estab- 
lished, it is strictly first order. More precisely for fc = 4, 
at a critical value p c — .6885(5) the system has a phase 
transition from a high p phase where there is a giant clus- 
ter with a finite density to a low-p phase where there is 
not. The density of the giant cluster is given by the fc-core 
density minus the density of the small clusters, that is ap- 
proximately Psmaii ~ 0.04 near the transition, therefore 
the critical properties of the giant cluster can be safely 
extracted from the total density in the percolating phase. 
The fc-core density exhibits a discontinuous transition, 
jumping from a p+ = 0.567(4) to p~ = 0.044(2). How- 
ever the density displays no singular behavior at the tran- 
sition, and the correlation length extracted from fc-core 
correlation function G(i,j) = {v%Vj) — (yj)(vj) (where V\ 
is 1 on the fc-core and otherwise [l(J]) remains finite, 
tc < 10 ; 

We started the numerical investigations considering 
hypercubic lattices with periodic boundary conditions 
(PBC). In fig. 0} we plot the density of the fc-core 
for a sample of size L = 320, corresponding to O(10 10 ) 
sites. The density of the fc-core has a discontinuous tran- 
sition at a p = .6869, where it jumps from a high-density 
percolating phase to a low-density non-percolating phase 
p w .05. The behavior appears to be consistent with a 
singular behavior at the transition but a careful study 
of the data in order to extract the critical p and the ex- 
ponent (3 shows some inconsistencies. Indeed the curve 
seems to be fit at best with the exponent j3 = 1/2 (the 
mean-field Bethc-lattice value) but with a value of the 
critical probability p ~ .6862 definitively lower than that 
at which the transition is actually observed. In order 
to assess if this behavior can be considered a finite-size 
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FIG. 1: fc-core density vs probability with k — 4 for a four- 
dimensional sample of size L = 320 and periodic boundary 
conditions. The arrow marks a discontinuous transition in 
the density that jumps to p » .05 at lower values of p. The 
data are fit in the region near the transition with the function 
p(p) ~ p(p c )+b(p— Pc) 1 ^ 2 but the actual transition probability 
(marked by the arrow) is above the p c estimated from the fit. 
Inset: rescaled plot of the data and the fit near the transition 



FIG. 2: Inverse of the correlation length £ vs. probability for 
the same sample of fig. Q The arrow marks the point of the 
actual transition. Although the behavior is consistent with a 
divergence at p ~ .6862, in all sample studied we never ob- 
served a correlation length at the actual transition exceeding 
£ ~ 8 — 10, i.e. large but much smaller than the sample size 
L = 320. The vertical line marks the value of the true critical 
probability p c , see text. 



effect we investigated the fc-core spatial correlation func- 
tion. In figure @ we plot the inverse correlation length 
£ as a function of p for the same sample, again the plot 
is apparently consistent with a divergence of £ at a value 
of p slightly lower than the value at which the transi- 
tion is actually observed (marked with an arrow in fig. 
|2J). In principle this could be a finite size effect but the 
problem is that the value of £ at the actual transition is 
large (£ ~ 8 — 10) but not comparable with the size of the 
system (L = 320). Actually we were not able to observe 
a correlation length bigger that ten lattice spacing inde- 
pendently of the size of the system (up to L — 320 ) and 
it seems highly unlikely that this divergence drives the 
transition. These features made us suspect that the the 
actual mechanism driving the transition is not the diver- 
gence of the correlation length but rather the nucleation 
of droplets of the low-density phase. Due to the periodic 
boundary conditions these nucleation centers are origi- 
nally absent in the system but they appear as spatial 
fluctuations of the density when the correlation length is 
sufficiently big leading to the transition. We tested this 
idea by putting some nucleation center {i.e. empty hy- 
percubes of size I) by hand in the sample and checked 
that this procedure shifts the transition at higher values 
of p although the size of the hypercubes (I = 20, 40) is 
small with respect to the size of the system and large 
with respect to the correlation length. According to this 
interpretation the percolating phase is unstable and can 
be observed only because nucleation centers in finite-size 
systems with periodic boundary conditions are extremely 
rare. 



In order to assess the validity of this interpretation and 
to determine whether there is a true percolation transi- 
tion at higher values of p we considered systems with 
completely empty boundaries. These boundary condi- 
tions guarantee that the system is completely isolated 
from the outside, as a consequence the density at finite 
size is a lower bound to the density in the thermody- 
namic limit. Furthermore it turned out that in this case 
numerical methods are rather safe for extrapolating the 
behavior of the infinite system at variance with the more 
delicate k > d case where finite-size effects are extremely 
large as mentioned above 9] . The possibility of choosing 
these boundary condition is a special feature of the case 
k < d because otherwise no site can resist culling if the 
boundaries of the hypercube are empty. 

In order to estimate the density at a given value of p 
we generated lattices of increasing size. In figure © we 
plot the results for the total density at various sizes 
L. The total density pi, is affected by the presence of 
the empty boundaries and is an increasing function of 
L. The monotonicity property allows to safely conclude 
from fig. |J3j that there is indeed a discontinuous tran- 
sition in the large L limit. We also measured the bulk 
density, i.e. the density of a smaller hypercube inside 
the sample whose boundaries are far enough from the 
surfaces. The bulk density provides a direct estimate of 
the density in thermodynamic limit and strengthens the 
conclusion that there is a discontinuous transition, see fig. 
13 The transition probability decreases with the sample 
size and tends to a critical value p c = .6885(5) that was 
estimated through extrapolation. We expect that the in- 
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FIG. 3: bulk density and total density vs. probability for 
different sample sizes, the average total density at given size 
is a lower bound to the bulk density because of the empty 
boundary conditions. 



terface between the low-density phase on the boundaries 
and the high-density phase in the bulk penetrates more 
and more in the sample for p — > p c , as a consequence in 
smaller systems the transition will occur at higher val- 
ues of p. The percolation value p c = .6885(5) is larger 
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FIG. 4: Density vs. probability for a sample with periodic 
boundary conditions and for the bulk of a sample with empty 
boundary condition with the same size L — 320. The data 
and the fit for periodic boundary conditions are the same 
of in fig. The vertical line marks the estimated critical 
probability p c — .6885(5). The two arrows mark the position 
of the actual transition probabilities, the transition with PBC 
is in the unstable region, see text. Inset: magnified view of 
the bulk density near the transition computed on internal 
hypercubes of various sizes (I — 100, 150, 200) inside a sample 
with L = 320. 



than the actual value of the transition in the case of pe- 



riodic boundary condition, e.g. p = .6869 for L = 320, 
in figure 0] we plot the bulk density and compare it with 
the density of the system in the case of periodic bound- 
ary conditions. For p > p c (marked by a vertical line in 
the plot) the two densities are equal and we clearly see 
that the density in the case of periodic boundary condi- 
tions is the analytic continuation of the percolating- phase 
density in the unstable region. We also verified that the 
correlation length of the bulk is the same of the systems 
with periodic boundary conditions at the corresponding 
values of p, see fig. By looking at figures ^ and |21 we 
immediately see that the density and correlation length 
are regular at the estimated value of the real transition 
probability p c = .6885. In particular using this value 
of p c we estimate £ c 2.5, while in the unstable phase 
we can observe £ up to ten lattice spacings. We note 
that due to time and memory constraints, few lattices of 
the largest size (L = 320) can be studied, however the 
sample-to-sample fluctuations are small enough to ensure 
a good estimate the density, for instance the bulk den- 
sity estimates computed from hypercubes of different size 
L = 100, 150, 200 do not show significant deviations, see 
inset of fig. [3 

The qualitative features of the finite size effect in the 
case of empty boundaries can be understood in the same 
way as in thermodynamical first-order transitions, e.g. a 
ferromagnetic Ising model in a small positive field with 
the spins on the boundaries forced to be negative. As we 
noted above, the critical transition probability is shifted 
to higher values of p and the true p c can be estimated 
through extrapolation (for instance for L = 320 the tran- 
sition is at p — .689, see [2J ■ In general the finite size 
corrections to the density should scale as 1/L, i.e. 
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Pl{p) = p{p) - jc^p) + O(j^) 



(1) 



and figure |3 suggests that the factor c\(p) diverges at p c . 
This factor is determined by the density profiles p(p, z) 
at distance z from an empty surface: 



Clip) 



(pip) ~ pip,z))dz . 



(2) 



Direct inspection of the profiles pip, z) shows that they 
are consistent with a divergence at the transition, con- 
sistently with the expectation that the transition is de- 
termined by the penetration deep inside the sample of 
the interface between the low-density and the percolat- 
ing phases. The precise nature of the divergence of C\ip) 
would require a more detailed analysis which goes beyond 
the scope of this work. In figure (JSJ we plot the behavior 
of the inverse of c\ and of the transition probability for 
different sample sizes from which p c was estimated. 

In conclusion the density of the fc-core in four dimen- 
sion with k = 4 exhibits a discontinuous transition at 
p c = .6885(5) from a high-density percolating phase 
to a low-density non-percolating phase. The transition 
loses its hybrid character with respect to the Bethe lat- 
tice case: the density is regular at the transition and 
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the correlation length is finite. These results are most 
clearly seen numerically considering samples with empty 
boundary conditions, instead in finite size systems with 
periodic boundary conditions it is possible to follow the 
percolating phase in the unstable region p < .6885(5). 
The behaviour of the unstable phase is consistent with a 
psettiio-transition at a lower probability p w .686. This 
pseudo-transition appears to have a hybrid character 
with pseudo-exponent (3 = 1/2 and diverging correla- 
tion length but it cannot be observed because the un- 
stable phase decays at p « .6869 through nucleation of 
low-density droplets determined by spatial fluctuations 
of the density. 



FIG. 5: probability vs. inverse of the 1/L correction ci for 
different sample sizes, the data are consistent with a diver- 
gence of the prefactors at the transition. Inset: probability 
at which the total density is p — .3 < p c as a function of the 
inverse of the size of the sample, it tends to p c as the size of 
the system increases. 



It would be very interesting to confirm these results 
through analytical methods. Although the 1/d expan- 
sion |10| gives no hint about the observed disappearance 
of the mixed transition, maybe different approaches to 
perturbation theory around the Bethe solution [l^. l2fj| 
could be able to see it, notwithstanding the possibility 
that non-perturbative effects should be taken into ac- 
count. 
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